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Introduction 

Fast space trips are important to intercept and rendezvous with an impacting asteroid or comet, 
particularly those not detected many years in advance. Fast trajectories can shorten space flight times and 
allow orbit modification efforts to begin earlier. The earlier the effort begins, the less change in velocity 
(AV) required to alter the object’s trajectory. However, shorter trip times require more propellant to 
provide enough thrust if the travel distance is fixed. This additional propellant mass can be a burden to 
the structural architecture of spacecraft. Thus, it is necessary to trade off between flight time and propel- 
lant mass. 

Background 

This report describes the approach and results of an end-to-end simulation to deflect a long-period 
comet (LPC) by using a rapid rendezvous spacecraft and laser ablation system. The laser energy required 
for providing sufficient deflection AV and an analysis of possible intercept/rendezvous spacecraft trajec- 
tories are studied in this analysis. These problems minimize a weighted sum of the flight time and 
required propellant by using an advanced propulsion system. The optimal thrust-vector history and 
propellant mass to use are found in order to transfer a spacecraft from the Earth to a targeted celestial 
object. One goal of this analysis is to formulate an optimization problem for intercept/rendezvous space- 
craft trajectories. 

One approach to alter the trajectory of the object in a highly controlled manner is to use pulsed laser 
ablative propulsion (ref. 1). A sufficiently intense laser pulse ablates the surface of a near-Earth object 
(NEO) by causing plasma blowoff. The momentum change from a single laser pulse is very small. 
However, the cumulative effect is very effective because the laser can interact with the object over long 
periods of time. The laser ablation technique can overcome the mass penalties associated with other 
nondisruptive approaches because no propellant is required to generate the AV (the material of the celes- 
tial object is the propellant source). Additionally, laser ablation is effective against a wide range of 
surface materials and does not require any landing or physical attachment to the object. For diverting 
distant asteroids and comets, the power and optical requirements of a laser ablation system on or near the 
Earth may be too extreme to contemplate in the next few decades. A hybrid solution would be for a 
spacecraft to carry a laser as a payload to a particular celestial body. The spacecraft would require an 
advanced propulsion system capable of rapid rendezvous with the object and an extremely powerful 
electrical generator, which is likely needed for the propulsion system as well. The spacecraft would 
station-keep with the object at a “small” standoff distance while the laser ablation is performed. 


7 Chapter nomenclature available in chapter notes, p. 217. 


148 


Trajectory Optimization Algorithm 

For intercept and rendezvous trajectories, optimization problems (ref. 2) in three dimensions are for- 
mulated to minimize flight time with moderate propellant mass. Many problems in the design of modem 
guidance and control systems require optimization of the trajectory, which minimizes (or maximizes) 
some performance criterion. Using the theory of the calculus of variations, the formulation of these prob- 
lems yields a two-point boundary-value problem (TPBVP). The resulting optimal trajectory also satisfies 
the physical constraints and the given differential equations. The open-loop optimal trajectory can be used 
as a reference trajectory for intercept or rendezvous with Earth-crossing asteroids or comets. To simplify 
the presentation and focus more on the inequality constraint, we first present the necessary conditions for 
an optimal control problem without the inequality constraint and then discuss the inequality constraint 
separately. A general optimal control problem can be stated as follows: 

Given the performance index, J, and radial velocity, u, 

rtf 

J(u) = 4>[x(t f ),tf ]+ L(x,«,t) dt (1) 

Jt o 

subject to the dynamic equations, x, and boundary conditions 

x= f(x,n,t), x(t 0 ) = x 0 , t 0 given (2) 

and with free final time tf, find the control history u{ t) to minimize J(u) with the prescribed terminal 
constraints 

O[x(t f ),t f ] = 0 (3) 

Here x(t) e R n are the state variables, u{ t) e R/ are the control components, and <X> e R k . The Hamiltonian 
function is defined with Lagrange multipliers A(t) e R n as 

HsL+A T f (4) 

The performance index in equation (1) is augmented and rewritten as 

J '= 4 ) [x(t f ),t f ] + v T 0[x(t f ),t f ] + [ tf (H- A T x) dt (5) 

Jt o 

where V is a constant multiplier vector of the dimension of the constraint <X>. The Minimum Principle 
requires that the optimal controls minimize the Hamiltonian function, H: 

«*(t) = min arg H(x*,A*,M,t) (6) 

ugQ 


where Q is the set of admissible controls and x*, A*, and it * are the extremal states, costates, and 
controls, respectively. The initial time and the initial states are kn own. The conditions to be satisfied to 
minimize J ' are found by taking the first variation of J ' and setting relations equal to zero. From this, 
the states, costates, and the Hamiltonian function satisfy the following conditions: 

x T = H a (7a) 
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(7b) 


A T =-H x 


A T (t f ) = (^ + v T O x ) 

(7c) 

1 A=t f 


H(t f ) = -4 t +v T O t ) 

(7d) 

v A=t f 


H ; ,=0 

(7e) 


Equation (7e) can be solved for the control, so that the control is removed from equations (7a) and (7b). 
There are n + k + 1 u nk nown values; n A’s; k v’s; and the final time. These values can be solved by using 
the k terminal constraints, equation (3), the n equations, equation (7c), and equaton (7d). 

The control inequality constraint is represented as 

M min — ^(1) — M max 0 ^ t ^ tf (8) 

Control variable inequality constraint is augmented to the cost function and additional necessary 
conditions are obtained as a result (refs. 3 and 4). The optimal trajectory is composed of two types of 
control — nominal control [ w min < u( t) < « max ] and boundary control [ u( t) = « min or « max ]. The nomi- 
nal control satisfies the same necessary conditions as the unconstrained problem. For boundary control, 
the inequality constraint becomes an equality constraint. Many classical problems in the calculus of 
variations treat constraints of this form very well. A new Hamiltonian function with control variable 
inequality constraint is redefined as 

H — H + p j(n — n max ) + p.2( M min ~ M ) (9) 


where 


Pl(w « max ) 0 ft i — 0 I 

M»min- M ) = ° H- 2 — 0 J 
The necessary conditions and controls for active constraints are 

i T = H A 
A = -H x 

h m =o 


( 10 ) 


( 11 ) 

( 12 ) 

(13) 


We use the control u( t) from the condition H„ = 0 when the control constraints are not active. For the 
problem with active inequality constraints, equations (11) and (12) are the equations of state and costate 
variables. For the Lagrange multiplier p;, i = 1,2 must necessarily satisfy that 

p i = 0 if the associated constraint is not active 
p ; > 0 if the associated constraint is active 
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and (i; can be obtained from solving H„ =0 for (j.j. Hence, when u(t) < u m i n , the control and the 
Lagrange multipliers are u(t) = u m i n , u | = 0, u 2 = H u . When u(t) > u max , the control and the Lagrange 
multipliers are u(t) = u max , tl 2 = 0, u | = -H u . 

Many numerical algorithms to solve optimal control problems have been developed. Indirect methods 
are theoretically based on the Minimum Principle, which characterizes the set of optimal states and con- 
trols in terms of the solution of a boundary value problem. One indirect method is the shooting method, 
which yields solutions of high precision. The shooting method is a second order method and hence is 
very sensitive to small changes of costate initial conditions. Shooting methods have the associated diffi- 
culties caused by instability of the initial value problem for the system of differential equations and by the 
requirement for good initial guesses for the iterative solutions of nonlinear problems. In this analysis, a 
shooting method is used to solve the comets or asteroids intercept or rendezvous trajectory problems. 

Problem Statement 

Propulsion 

Many future propulsion systems have been proposed and analyzed. One potential propulsion approach 
that has been examined for a Comet/Asteroid Protection System (CAPS) deflection capability is the 
Variable Specific Impulse Magnetoplasma Rocket (VASIMR). VASIMR is a high power magnetoplasma 
rocket that gives continuous and variable thrust at constant power (ref. 5). Hydrogen plasma is heated by 
radio frequency (RF) power to increase exhaust velocity up to 300 km/s. The power output of the engine 
is kept constant, thus thrust and specific impulse, I sp , are inversely related. Thrust is increased propor- 
tional to the power level. The engine can optimize propellant usage and deliver a maximum payload in 
minimum time by varying thrust and I sp (ref. 6). Therefore, VASIMR can yield the fastest possible trip 
time with a given amount of propellant by using constant power throttling (CPT). AlO-kW space demon- 
strator experiment has been completed, and a VASIMR engine with 200-MW power could be available 
around the year 2050. The specific impulse range of the engine would be 3000 s to 30000 s, and the 
corresponding thrust range would be approximately 13600 N to 1360 N (assuming 100 percent power 
efficiency of 200 MW). To calculate acceleration, a, and spacecraft mass flow rate, the following rela- 
tionships are used. The thrust, T = |m|v e , and exhaust velocity, v e = I sp g 0 , are described by specific 
impulse, I sp , and the acceleration due to gravity at the Earth’s sea level, g 0 . Mass flow rate, m, is itself 
negative value. 


1 I . I 9 

The power, p, required to expel mass at the mass flow rate, m, is ep = — |m| v e . 8 is the efficiency of 


the propulsion system. Thus, we know 

1 ~ ^ 2ep 2ep 

2 v e ispgo 

Using equation (14), acceleration due to VASIMR can be derived: 


(14) 


T 2ep I 
a = — = - 


m mg G I 


sp 


(15) 
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where m is spacecraft mass at any time. Finally, the mass flow rate is calculated as 


-T 2ep 1 

V e " go Isp 


(16) 


Equations of Motion 

The spacecraft is considered to fly in three-dimensional interplanetary space. The three-degree-of- 
freedom equations of motion are the following: 


r = u 

(17a) 

• y2 W 2 ■ n 1 

u= — H + a sin a cos p — y 

r r r 2 

(17b) 

uv vw sin (b n 

v = 1 1- a cos a cos (3 

r r cos (f> 

(17c) 

uw v 2 sin (b . _ 

w = + a sm p 

r r cos <() 

(17d) 

9= V 
r cos ([) 

(17e) 

; w 

c^ = — 
r 

(17f) 

2ep 

m_ 2 t 2 

(17g) 


where r is the radial distance from the Sun to spacecraft, u is the radial velocity, v is the tangential 
velocity, w is the normal velocity, 9 is the angle measured from the x-axis (defined as vernal equinox) in 
the x-y plane, (]) is the angle measured from x-y plane, m is the mass of spacecraft, a is acceleration of 
spacecraft, p is the power of spacecraft, g is the gravitational parameter at Earth’s sea level, and I sp is the 
specific impulse of spacecraft engine. The control variables are thrust direction angle in plane (a), thrust 
direction angle of out-of plane ((3), and the specific impulse (I sp ). 

The Hamiltonian function is 


1 1 = w t + A r u+A u 


(22 , 1 
V W . Q 1 

— -I + a sin a cos p — - 


\ 


+ A, 


7 


^ uv vw sin (b n ^ 

1 1- a cos a cos p 

r r cos (|) 


+A V 


f 9 

uw v sin ()) 

r r cos (|) 


+ a sin (3 


+ Ac 


r cos {p T 


w 

+ Aa— -A 


2ep 

m 7 T~ 
go Isp 


( 18 ) 
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The costate equations are 


A r = A u 


f 2 2 

V + W 

2 ' 

+ A v 

vw sin 4> 

h 

uv 

+ A W 

^ 2 

v sin (f) 

\ 

UV 

< 2 

r 3 J 

v r 2 cos 4> 

r 2 J 

r 2 cos <j> 

r 2 ) 


+ A 


0 2 . 
r cos q) 


» w 
r 


(19a) 


V w 

A u = — A r + A v — h A w — 

r r 


A v =-A u — + A V 

r 


( u w sin <|) ^ 

r r cos ()) 


. 2v sin (b . 1 

+ A W ~ “Ae X 

r cos (p r cos tp 


(19b) 

(19c) 


« » 2w . v sin (() . u . 1 

- _A -u A v T + A 'd> - 

r r cos q) r Y r 


A , 0 -0 


Aa - —A, 


VW 


-+ A,,; 


--Ac 


v sin (|) 


r cos (|) r cos (j) r cos § 


(19d) 

(19e) 

(19f) 


2ep 2ep 2ep 

A m = A„ — ^ sin a cos p + A v — ~ cos a cos p + A 


m "u 2 , 
m gl. 


sp 


v 2 , 
m gl 


sp 


w 2 , 

m gl 


-sin (3 


sp 


(19g) 


Performance Index 


Optimal control theory is concerned with finding the control history to optimize a measure of the per- 
formance index of the following general form: 


f tf 

J (u) = -m t - + w t dt 

•’tn 


( 20 ) 


where mf represents the final mass of spacecraft, and w t is weight for flight time (w t is set as 10 in this 
analysis). For the problem at hand, it is required to find the optimal trajectory that maximizes the final 
mass of spacecraft and minimizes the flight time. 


Initial Conditions and Terminal Conditions 


The spacecraft departs from the Earth with the following initial conditions at t = 0 s: 

r(to) = 1 au, u(to) = 0, v(to) = Earth’s velocity, w(to) = 0, 

9(to) = obtained from spacecraft departure time before collision, 

9(to) = 9°, m(to) = free (unknown), tf = free (unknown) 

The initial mass of spacecraft is a free parameter to include propellant mass. Thus, initially the costate of 
mass is set as A m (to) = -1. The other costate values [A r (to), A u (to), A v (to), A w (to), Ae(to), A^to)] at to are 
u nkn own. 
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Final boundary conditions are specified to satisfy the position and velocity of asteroids/comets. The 
spacecraft must intercept or rendezvous with the targeted asteroid/comet with specified orbit. To achieve 
the desired trajectory, final conditions should be satisfied. These are the positions of spacecraft for inter- 
cept trajectory and the positions and velocities of spacecraft for rendezvous trajectory. Hence, for inter- 
cept trajectory, the terminal state conditions are 


0[x(t f ),t f ]s< 


r (tf ) r target (tf ) 
9(tf ) — ^target (*f) 
) — ^target (tf ) 


( = 0 


m (tf)-mdry 


( 21 ) 


Rendezvous trajectory has the following terminal state conditions: 


0[x(t f ),t f ] = 


r (tf ) — Target (tf ) 

U(tf) _ Ujgj-ggf (tf ) 

v (lf ) — v target (tf ) 

w(tf)- w target (t f ) 

0(lf ) — ^target (tf) 
^(tf ) — ^target (tf ) 
m (t f )-m dry 


1 = 0 


( 22 ) 


where subscript “target” denotes the state of targeted celestial object, and m dry is the spacecraft dry mass. 
From transversality conditions, we obtain the following costate terminal conditions for intercept 
trajectory: 


Mtf) = vi 

(23a) 

Mtf) = Mt f ) = x w (tf) = o 

(23b) 

^e(tf) = v 2 

(23c) 

^(|>(tf) = v 3 

(23d) 

^m(tf) = v 4 

(23e) 


where V; is the Lagrange multiplier. From transversality conditions, we obtain the following costate 
terminal conditions for rendezvous trajectory: 


^r(tf) = Vi 

(24a) 

^'u( t f) = v 2 

(24b) 

^'v( t f ) = v 3 

(24c) 
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A, w (tf) = V 4 

(24d) 

^e(tf) = v 5 

(24e) 

^(|)( t f) = v 6 

(241) 

^m( t f) = v 7 

(24g) 


Furthermore, from equation (7d) the following condition is also satisfied at tf: 

H(t f ) = 0 (25) 

Equation (25) becomes another boundary equation at the final time. For the intercept problem, there 
are 14 differential equations describing the states and costates, with 15 unknowns [m(to), A, r (to), A u (to), 
^v(to), A-w(to), ^e(to), A-q(to), tf, Vi, v 2 , v 3 , v 4 , u(tf), v(tf), and w(t f )] and 15 boundary conditions [r(t 0 ), 
u(to), v(to), w(t 0 ), 0(t o ), <|>(to), A-m(to), eqs. (21), A, u (t f ), A, v (t f ), A. w (t f ), and eq. (25)]. For the rendezvous 
problem, there are 14 differential equations describing the states and costates, with 15 unknowns [m(to), 
^r(to), A, u (t 0 ), A, v (to), A, w (to), ^e(to), ^q(to), tf, Vi, v 2 , v 3 , v 4 , v 5 , v 6 , and v 7 ] and 15 boundary conditions 
[r(to), u(to), v(to), w(to), 0(to), <f>(to), X m (to), eqs. (22), and eq. (25)]. Thus, the two-point boundary prob- 
lem can be completely solved with these boundary conditions. 

Controls 

The control variables are thrust direction angle in-plane (a), thrust direction angle of out-of-plane ((3), 
and specific impulse (I sp ). A second-order necessary condition, the Legendre condition, states that the 
second derivative of the Hamiltonian, with respect to the controls, must be greater than or equal to zero 
for the performance index to be at a minimum. Thus, H aa , Hpp must be greater than or equal to zero for 
the performance index to be at a minimum. The first derivative of H with respect to a and convexity 
condition yields a control variable of a as follows: 


sin a = 




cos a = 


k u + ^ V 


(26a) 


(26b) 


The first derivative of H with respect to (3 and convexity condition yields a control variable of (3 as 
follows: 


sin (3 = 


-A. v 


■\J k u + k y + k 


2 

w 


cos (3 = 


k u + ^ V 


k~„ + k v + k~. 


(27a) 


(27b) 
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The first derivative of H with respect to I sp and convexity condition yields a control variable of I sp as 
follows: 


-2rnX m 

Sp U 2 , ,2 , ,2 


§oA/^u + ^v 


1 <1 <1 

A spmin — A sp — A spmax 


(28) 

(29) 


If the I sp constraint is not active and A, m < 0, the solution would be optimal. If the unconstraint control 
Isp ^ Isp min; then I sp — I sp min; whereas I sp — I sp max if the unconstraint control I sp > I sp max- for ^ 0; 
the optimal solution is going to be I sp = I sp m i n . 

Interplanetary optimal trajectories have been computed to maximize the final mass of the vehicle and 
to minimize the flight time. The optimal problem is a free final time problem to find the three controls 
satisfying the state and costate equations. 

End-to-End Simulation 


Here we consider a fictitious impacting LPC whose orbital parameters are given by aphelial distance 
r a = 100 au and perihelion distance r p = 0.7 au with inclination i = 50°. These orbital parameters yield 
semimajor axis a = 50.35 au, eccentricity e = 0.98609732, and orbital period of 357.27 years. It is 
assumed that the LPC has its density as p = 1000 kg/m 3 . The minimum required impulses for deflecting 
an impactor by 3 Earth radii are solved in this analysis, and the calculation is always performed to move 
the LPC’s trajectory from crossing the Earth’s orbit at the Earth’s center. The solutions represent impulse 
vectors that can be described by the magnitude of the minimum impulse and the optimal impulse angle. 
The gravitation effects of Earth are considered by using a three-dimensional optimization problem to 
calculate the impulse vectors. Figure 1 includes the minimum AV with respect to the impulse time that is 
defined as time before collision when AV is applied. Figure 2 includes the optimal impulse angle with 
respect to the impulse time. The impulse angle is described as being in an asteroid’s/comet’s orbital plane 
and is defined as the angle from the asteroid’s/comet’s original velocity vector to the impulse vector 
toward the Sun-asteroid/comet line. The dotted line in these figures explains the preperihelial collision 
case (a collision occurs before an asteroid/comet passes its perihelion), whereas the solid line in these 
figures explains the postperihelial collision case (a collision occurs after an asteroid/comet passes its 
perihelion). Figure 3 shows an estimate of the typical energy required for laser ablation to deflect the 
1-km comet by 3 Earth radii, when preperihelial collision with Earth is considered. The required laser 
energy is also a function of the object’s density and the required AV, which varies depending on the 
object’s orbit and when the deflection occurs. It is easy to estimate the required laser energy for any size 
comet (or asteroid) because the minimum AV is linearly proportional to the cube of its diameter. Figure 4 
explains a detail of figure 3 for an LPC, using fixed impulse times, and describes the required energy for 
laser ablation to deflect the given size of a comet by 3 Earth radii. The estimated energy is calculated 
assuming that the cumulative energy generated by the laser is applied as an equivalent impulsive AV at 
some time before collision. Because the laser ablation occurs over a significant period of time, the laser 
interaction must be complete prior to the time specified for each curve in order to assure that the deflec- 
tion could be accomplished. Figure 5 shows a preliminary estimate of the achieved energy for a given 
laser power and operation period. Figures 4 and 5 can be used to estimate the nominal laser power 
required for a deflection mission. To illustrate a preperihelial case, figure 4 shows that approximately 
5 x 10 4 GJ of energy is required to deflect a 0.2-km comet by 3 Earth radii if applied 1 year before 
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Impulse angle, deg Minimum AV, cm/s 



Figure 1. Minimum AV to deflect 1-km LPC by 3 Earth radii. 



Figure 2. Optimal impulse angle with respect to impulse time. 
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Figure 3. Required laser energy to deflect 1-km LPC by 3 Earth radii. 



Figure 4. Estimated energy required for laser ablation versus diameter for LPCs (p = 1000 kg/m 3 ). 
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Figure 5. Estimate of achieved energy for given laser power and operation period. 

collision. Figure 5 shows that a 100-MW laser (10000-kJ pulse and 10-Hz laser repetition frequency) 
would be required to operate continuously for approximately 5 days in order to achieve this cumulative 
energy. Figure 6 illustrates the trajectories and thrust vector for a 100-metric ton (t) spacecraft (including 
10-t payload) to intercept or rendezvous with this particular LPC for a departure time of 7 months before 
a preperihelial collision with Earth. 

The propulsion system is assumed to operate at 90-percent power efficiency. Figures 7 and 8 show 
flight time and required propellant for intercept and rendezvous trajectories for various departure times 
when a 100-t spacecraft with 1 GW of power is assumed. For this specific LPC, there is a peak at 
11 months departure time because the spacecraft must fly in the reverse direction with respect to Earth’s 
orbital velocity. Because there are local minima and maxima in the propellant required and flight time, as 
shown in the figures, it can be concluded that the values are dependent upon the orbital geometry relation- 
ship as well as distance between the Earth and the comet. For a given departure time, the rendezvous 
trajectory requires more propellant and longer flight time than the intercept trajectory. This is because the 
terminal velocity of rendezvous spacecraft must be matched with the target’s velocity, which is not 
required for the intercept trajectory. Asteroids and comets with different orbital elements will have 
different flight times and propellant requirements. Even for the same celestial object, a postperihelial 
impact would have different results from those of a preperihelial impact. 

Once the spacecraft has rendezvoused with the LPC, a laser ablation system makes use of the same 
electrical power system that the propulsion system uses for the orbital transfer. For example, if the laser 
ablation operation can be completed 12 months before collision, approximately 3 x 10 6 GJ of energy 
would be required for a 0.8-km LPC to be deflected by 3 Earth radii. If we choose a 500-MW laser 
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Flight time to comet, months 



Figure 6. Example long-period comet intercept and rendezvous trajectories. 
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Figure 8. Propellant used for each departure time. 


system (assuming 1 GW of supplied power and 50-percent laser system efficiency), it would take about 
90 days of continuous operation to provide this amount of energy. From figure 7, the rendezvous space- 
craft with a laser ablation system should depart from Earth approximately 19 months before the collision, 
with a trip time of about 4 months. For the mission, approximately 340 t of propellant would be required, 
assuming 100-t total dry mass of the spacecraft and payload. Depending on the payload mass capability 
of the propulsion system, multiple deflection devices could be delivered to the impactor, including a 
fallback option to the laser ablation system. In case the rendezvous deflection mission was unsuccessful, 
a similar spacecraft with a different payload (e.g., nuclear explosive device) could be sent to deflect the 
target using an intercept trajectory. If we assume that the intercept spacecraft departs from Earth 
9 months before the collision with a 10-t payload, the spacecraft can arrive at the target approximately 
7 months before impact and may require less than 230 t of propellant (as shown in figs. 7 and 8). 

A multimegawatt nuclear electric propulsion (NEP) system using a YASIMR engine is currently esti- 
mated to have a maximum overall specific mass of 1.0 kg/kW (ref. 7). For a 1-GW system, this would 
result in a total spacecraft dry mass of 1000 mt (neglecting payload mass). This is ten times the total 
spacecraft mass assumed for this analysis, and the mass of a future laser ablation payload is presently not 
well understood. The assumed power efficiency of 90 percent is also optimistic. More capable power 
generators (gigawatt class) with lower specific masses could provide the power needed to reduce trip 
times and provide more powerful lasers. Deflecting an impactor by only 3 Earth radii would be sufficient 
for a deflection effort. More powerful lasers would be capable of providing a greater miss distance, and 
thus more margin for uncertainty in the object’s orbit. NEOs with greater densities would also require a 
more capable laser ablation system. Longer warning times would reduce the requirements on the orbit 
modification system but would make the CAPS detection system more challenging to implement. For 
the 0.8-km LPC assumed in this analysis, the detection system would need to determine the comet’s 
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trajectory at least 19 months prior to its preperihelial collision with the Earth. For impactors with 
extremely short warning times, an intercept trajectory may be the only feasible scenario for diverting the 
object. 

Technological advances that can significantly reduce the specific mass of the rendezvous spacecraft 
and laser payload may permit this type of deflection approach to become a reality. A tiered planetary 
defense approach using rapid rendezvous and intercept spacecraft could provide a feasible scenario to 
protect the Earth from an impacting LPC, as well as other classes of impacting NEOs. A rendezvous 
spacecraft with a laser ablation payload could also provide a capable and robust orbit modification 
approach for altering an NEO’s orbit for resource utilization. 

Concluding Remarks 

This report presents intercept/rendezvous trajectories for an advanced spacecraft that is designed to 
deliver laser ablation energy to an Earth-crossing long-period comet. The trajectory optimization problem 
is solved using the shooting method, which yields highly accurate solutions. The open-loop optimal 
solutions can be used as reference spacecraft trajectory for the deflection problem. The end-to-end simu- 
lation in this report demonstrates a conceptual approach to altering the orbit of an Earth impacting long- 
period comet, particularly one which represents an immediate threat. 
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